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Abstract. The efficiency of star formation governs many observable properties of the cosmologi- 
cal galaxy population, yet many current models of galaxy formation largely ignore the important 
physics of star formation and the interstellar medium (ISM). Using hydrodynamical simulations 
of disk galaxies that include a treatment of the molecular ISM and star formation in molecular 
clouds (Robertson & Kravtsov 2008), we study the influence of star formation efficiency and 
molecular hydrogen abundance on the properties of high-redshift galaxy populations. In this 
work, we focus on a model of low-mass, star forming galaxies at 1 < z < 2 that may host long 
duration gamma-ray bursts (GRBs). Observations of GRB hosts have revealed a population of 
faint systems with star formation properties that often differ from Lyman-break galaxies (LBGs) 
and more luminous high-redshift field galaxies. Observed GRB sightlines are deficient in molec- 
ular hydrogen, but it is unclear to what degree this deficiency owes to intrinsic properties of the 
galaxy or the impact the GRB has on its environment. We find that hydrodynamical simulations 
of low-stellar mass systems at high-redshifts can reproduce the observed star formation rates 
and efficiencies of GRB host galaxies at redshifts 1 < z < 2. We show that the compact structure 
of low-mass high-redshift GRB hosts may lead to a molecular ISM fraction of a few tenths, well 
above that observed in individual GRB sightlines. However, the star formation rates of observed 
GRB host galaxies imply molecular gas masses of 10 s — 10 9 Mq similar to those produced in 
the simulations, and may therefore imply fairly large average H2 fractions in their ISM. 
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1. Introduction 

To improve the physical description of star formation in hydrodynamical simulations 
of galaxies, Robertson & Kravtsov (2008) implemented a new model for the ISM that 
includes low-temperature (T < 10 4 K) cooling, directly ties the star formation rate to 
the molecular gas density, and accounts for the destruction of molecular hydrogen by an 
interstellar radiation field (ISRF) from young stars. They used simulations to study the 
relation between star formation and the ISM in galaxies and demonstrated that, for the 
first time, their new model simultaneously reproduces the molecular gas and total gas 
Kennicutt-Schmidt (KS) relations, the connection between star formation and disk rota- 
tion, and the relation between interstellar pressure and the fraction of gas in molecular 
form (e.g. Wong & Blitz 2002, Blitz & Rosolowsky 2006). The capability of this model to 
reproduce both the star formation efficiency and molecular abundance of nearby systems 
makes it useful for simulating low-mass galaxies that have suppressed H2 abundances 
(and whose star formation rates would be overestimated in common treatments of star 
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formation based on the KS relation) and high-rcdshift galaxies whose structural prop- 
erties may vary substantially from local systems (and may therefore not have the same 
KS relation normalization). The model should be especially useful for studying low-mass 
galaxies at high-redshift, such as long duration gamma-ray burst (GRB) host galaxies at 
1 < z < 2, which is the focus of this work. 

The highly-energetic phenomena known as GRBs were discovered over forty years ago 
(Klcbcsadcl ct al. 1973), but their extragalactic origin was confirmed only in the last 
decade (e.g., Metzger et al. 1997). Since then, the properties of the cosmological popu- 
lation of galaxies that host GRBs have been increasingly well-studied (e.g., Bloom et al. 
2002, Le Floc'h et al. 2006, Prochaska et al. 2006, Berger et al. 2007a,b). Recently, interest 
in long duration GRB galaxy hosts as possible tracers of the global star formation history 
of the universe has motivated systematic studies of their star formation efficiencies and 
stellar masses (Castro Ceron et al. 2008, Savaglio et al. 2008). These studies have found 
that high-redshift GRB hosts have small stellar masses (logM* ~ 9.3) and moderate 
star formation rates (SFR <~ 2.5 M yr -1 ). Compared with other high-redshift galaxy 
populations, GRB hosts tend to have lower star formation rates at fixed stellar mass 
compared with Lyman-break galaxies and lower stellar masses at fixed star formation 
rate compared with field galaxies (for details, see Savaglio et al. 2008). 

Spectroscopic studies of GRB sightlines have provided additional information about 
the post-explosion character of the host galaxy ISM. Tumlinson et al. (2007) failed to 
detect H2 in five GRB sightlines and suggested that low metallicity and large far ultra- 
violet ISRF strengths (10 — 100 x the Milky Way value) were responsible for destroying 
molecular hydrogen in GRB hosts. They interpreted the lack of vibrationally excited H 2 
lines as evidence against the GRB destroying its parent molecular cloud, but noted var- 
ious caveats to this conclusion such as the parent cloud size or cloud photodissociation 
before to the GRB. Whalen ct al. (2008) used one-dimensional radiative hydrodynamical 
calculations to show that GRBs can ionize nearby neutral hydrogen, but suggested that 
an additional ISRF is necessary to remove molecular hydrogen from the nearby ISM. 
Prochaska et al. (2008) studied NV absorption in GRB sightlines, and argued that if 
nitrogen ionization by GRB afterglows leads to NV absorption then the observations 
support a scenario where dense, molecular cloud-like environments serve as the sites of 
GRBs. 

Given the increasingly detailed studies of GRB hosts, their interesting ISM and star 
formation properties, and their low stellar masses, a theoretical study of GRB host galaxy 
analogues using hydrodynamical simulations that include a treatment of the molecular 
ISM is warranted. Below, we present simulations of a model GRB host galaxy that 
include a prescription for the molecular ISM and star formation in molecular clouds 
(Robertson & Kravtsov 2008). We use the simulations to examine the star formation 
efficiency and molecular hydrogen content of galaxies with structural properties similar 
to those expected for low-mass galaxies at 1 < z < 2. Below, we discuss our methodology 
and present some initial results. 



2. Methodology 

To study the properties of long duration GRB host galaxies, we simulate a numerical 
model of an isolated galaxy using a version of the N-body/Smoothcd Particle Hydrody- 
namics code GADGET (Springcl ct al. 2001, Springcl 2005b) that incorporates a model 
for the molecular ISM (Robertson & Kravtsov 2008). For details regarding the numerical 
galaxy models, simulation methodology, and ISM model, we refer the reader to Springcl 
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Figure 1. Simulated long duration Gamma-Ray Burst (GRB) host galaxy analogue at z ~ 1.5. 
Shown is the gas surface density of the GRB host (image intensity), color coded by the me- 
dian interstellar medium temperature (purple regions have T < 10 3 K, while blue regions have 
T > 10 4 K). The simulated galaxy has a stellar mass logAf* w 9.3 and a star formation rate 
SFR w 1.2M0jrr — 1 , similar to high-redshift GRB host galaxies (e.g., Castro Ceron et al. 2008, 
Savaglio et al. 2008). The simulations include the Robertson & Kravtsov (2008) model of the 
molecular ISM, enabling a study of the connection between star formation rate, galaxy proper- 
ties, and H2 abundance in GRB hosts. 

et al. (2005a), Robertson et al. (2006a, b), and Robertson & Kravtsov (2008), but a brief 
summary follows. 

The numerical galaxy model is designed to approximate the properties of 1 < z < 2 
GRB host galaxies as determined by Savaglio et al. 2008. The stellar disk mass of the 
system is set to logM* = 9.3, with a gas fraction of / gas = 0.5 (appropriate for high- 
redshift, see Erb et al. 2006), which implies a total virial mass of logM V i r = 10.9 for 
a typical disk baryon fraction of /b = 0.05. The virial radius is set appropriately for 
a halo with virial mass M v - U at z ~ 2. The exponential disk scale length was fixed 
according to the Mo et al. (1998) formalism, including the adjustment for an effective 
Navarro et al. (1996) dark matter halo concentration of cnfw = 6 (also appropriate for 
the chosen virial mass and redshift, see Bullock et al. 2001) and a spin of A = 0.05. 
The density field of the dark matter halo follows the Hernquist (1990) profile, while the 
velocity fields of the dark matter halo and the exponential stellar disk are set using the 
Hernquist (1990) distribution function and the epicyclical approximation, respectively. 
The numerical realizations of the stellar disk, gaseous disk, and dark matter halo are 
initialized with N disk ^ = 4 x 10 5 , -/V d i S k,gas = 4 x 10 5 , and iV D M = 4x 10 5 particles, and 
are evolved with a gravitational softening of e = 70 pc. The simulation is calculated for 
a duration of t ~ 1 Gyr, or about the time between redshift z ~ 2 and z ~ 1.5. 

The simulation includes a treatment of the physics of the ISM and star formation 
following the model presented by Robertson & Kravtsov (2008), and interested readers 
should examine that work for details. The photoionization code CLOUDY (Ferland et 
al. 1998) is used to tabulate the cooling rate, heating rate, molecular abundance, and 
related properties of g 0jS CIS Ob function of density, temperature, metallicity, and local 
interstellar radiation field (ISRF) strength. The star formation rate is calculated by 
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Figure 2. Kennicutt-Schmidt relation for a simulated GRB host galaxy at z ~ 1.5. Shown is 
the star formation rate surface density Esfr as a function of total gas surface density £ gas , 
measured in annuli (blue dots). The average Kennicutt-Schmidt relation of the GRB host has 
a steeper power-law index (a ~ 3.0, green line) than the disk-averaged relation measured by 
Kennicutt (1998; a ~ 1.4, dashed line), owing to the suppression of H2 in the galaxy exterior 
by the interstellar radiation field and the low ISM metallicity (for a detailed discussion, see 
Robertson & Kravtsov 2008). 

converting the molecular gas density to stars on a timescale that scales with the local 
dynamical time, with an efficiency set to match the star formation efficiency per free fall 
time in local molecular clouds (e.g., Krumholz & McKee 2005, Krumholz & Tan 2007). 
The local ISRF spectral shape is fixed to the local Milky Way ISRF inferred by Mathis 
et al. (1983), but the ISRF strength scales with the local star formation rate density 
(i.e., young, massive stars supply the local ultraviolet radiation field). The abundance of 
molecular gas tracked using CLOUDY includes the photodissociative and heating effects 
of this ISRF, and thereby includes a coarse accounting of the regulatory impact of the 
ISRF on star formation in molecular clouds. 

3. Results 

Figure [1] shows the gaseous structure of the GRB host galaxy model at z ~ 1.5 (after 
900 Myr of evolution). The figure shows the gas surface density of the system (image 
intensity) and the median temperature of the local ISM (purple regions have temperatures 
T < 10 3 K, while blue regions have T > 10 4 K). The system has a rotational velocity of 
v mt « 100 km s _1 and a disk scale length of Rd ks 1.5 kpc. During the simulation 
the star formation rate of the system varies in the range SFR w 0.5 — 2.5 Mq yr _1 , 
while the specific star formation rate is SFR/M* rs 0.17 — 1.1 Gyr -1 . These properties 
are consistent with the properties of high-redshift GRB host galaxies determined by 
Savaglio et al. (2008), who find star formation rates of SFR ~ 0.1 — 10 and specific star 
formation rates of M*/SFR ~ 0.1 - 10 Gyr. 
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The compactness of the system leads to a dense ISM and a considerable molecular 
gas fraction. The global, mass-weighted molecular abundance declines from /h2 ~ 0.5 at 
z ~ 2 to /h2 ~ 0.3 at z ~ 1.5. As a function of radius, the molecular fraction declines from 
/h2 ~ 1 near the center of the galaxy to /h2 ~ 0.1 — 0.3 beyond a disk scale radius. The 
typical star formation rate-weighted radius of the system is rgpR ~ 0.8 kpc, where the 
ISM molecular fraction is /h2 ~ 0.6. Hence, if observed high-redshift GRB host galaxies 
are similar in nature to this simulated system, GRBs will likely occur in molecular-rich 
regions. While these molecular abundances are consistent with the spectroscopic studies 
by Prochaska et al. (2008), they are well above the observed H2 abundance along GRB 
sightlines (e.g., Tumlinson et al. 2007). In this model, the compact and dense structure 
of the high-redshift GRB host prevents the diffuse ISRF from suppress the H2 to levels 
observed in GRB sightlines. If an ISRF is responsible for suppressing H2 to observed 
levels in GRB sightlines, it may be generated by discrete point sources nearby the GRB 
in a manner not captured by the diffuse ISRF included in these simulations. 

We note that in order to supply the star formation efficiency for GRB hosts determined 
by Savaglio et al. (2008), GRB hosts may need to be fairly molecule rich if their structure 
is similar to the simulated high-redshift galaxy analogues presented here. Figure [5] shows 
the total gas Kennicutt-Schmidt (KS) relation for the GRB host galaxy analogue, mea- 
sured in annuli with a width of Ar — 100 pc. Plotted is the star formation rate density 
Esfr as a function of the total gas surface density E gas (blue points), compared with 
the mean disk-averaged trend determined by Kennicutt (1998, dashed-linc) . The central 
concentration of molecular gas causes the total gas KS relation of the simulated GRB 
host galaxy analogue to be steeper than the disk-averaged relation. In order to supply 
the observed star formation rate of SFR ~ 1 — 10 M©yr -1 , as this simulated galaxy does, 
the typical consumption timescales of ~ 100 Myr for molecular gas imply a reservoir of 
roughly M m - 0.1 - 1 x 10 9 M (the simulated system has M H2 ~ 2 - 7 x 10 8 M© 
during its evolution). Since observed GRB hosts have stellar masses of only logM* ~ 9.3 
(Savaglio et al. 2008), the inferred molecular fraction of the ISM should be large even for 
very gas rich systems. 

Overall, we find that under standard assumptions about the mass and redshift scalings 
of galaxy structure hydrodynamical simulations of disk galaxies with stellar masses of 
logAf* ~ 9 that utilize a model for the molecular ISM and star formation in molecular 
clouds (Robertson & Kravtsov 2008) can reproduce the observed star formation rates 
and efficiencies of GRB hosts (e.g., Castro Ceron et al. 2008, Savaglio et al. 2008). 
The star formation in both observed GRB hosts and the simulated GRB host analogue 
presented here is efficient for their low stellar masses, and to supply the observed range 
of star formation rates (SFR ~ 0.1 — 10 M© yr _1 ) the molecular gas content such 
systems may need to be considerable (/h2 > 0.1). While this result is consistent with 
observations that suggest GRBs occur in dense, potentially molecular-rich regions of the 
ISM (e.g. Prochaska et al. 2008), more work is needed to reconcile such results with 
the low molecular abundance observed in GRB sightlines (e.g. Tumlinson et al. 2007) if 
GRBs cannot efficiently destroy H 2 in the ISM (Whalen et al. 2008). 
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